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Using first-principles calculations, we show that the conduction and valence band energies and 
their deformation potentials exhibit a non-negligible compositional bowing in strained ternary semi¬ 
conductor alloys such as InGaAs. The electronic structure of these compounds has been calculated 
within the framework of local density approximation and hybrid functional approach for large cubic 
supercells and special quasi-random structures, which represent two kinds of model structures for 
random alloys. We find that the predicted bowing effect for the band energy deformation poten¬ 
tials is rather insensitive to the choice of the functional and alloy structural model. The direction 
of bowing is determined by In cations that give a stronger contribution to the formation of the 
In^Gai-^As valence band states with x > 0.5, compared to Ga cations. 


III-V semiconductor compounds are considered for the 
integration with Si-based microelectronics to take advan¬ 
tage of their high charge carrier mobility^ This has been 
driving basic and applied research on various aspects of 
the III-V semiconductor materials as well as on the device 
designA— One field of extensive research is to understand 
how the electronic properties of the III-V compounds are 
altered by stress that is externally applied to tailor char¬ 
acteristics of semiconductor nanodevices^^^ 

Though III-V compounds such as InGaAs are impor¬ 
tant material systems, measuring their electronic prop¬ 
erties is a formidable challenge, partly due to their nar¬ 
row gap. As a consequence, rather scattered experimen¬ 
tal data have been reported so far, especially for the X 
and L conduction band valleys^ suggesting that reli¬ 
able theoretical studies are needed. Knowing deforma¬ 
tion potentials of the band gap, valence band splitting 
and higher energies of conduction band minima is of par¬ 
ticular interest for opto- and nanoelectronics^ So far, 
the bulk of the theoretical work has been done within 
the framework of the empirical pseudopotential method 
combined with virtual crystal approximation (VCA)i£ 
These methods require adjustable parameters fitted to 
the existing experimental data. Nowadays, the state- 
of-the-art first-principles methods allow for large-scale 
atomistic, parameter-free calculations of the properties 
of disordered materials^ as well as parameter-free VGA 
calculations^ A recent first-principles study has pre¬ 
dicted a strong compositional bowing of deformation po¬ 
tentials in group Ill-nitride semiconductor compounds^ 
not supporting an assumption of a negligible bowing* 12 
There has been no systematic atomistic study of the com¬ 
positional bowing effect in the strained InGaAs ternary 
alloy, which is a channel material of high relevance for 
modern types of field-effect transistors A— 

In this Letter, we have done first-principles calcu¬ 
lations of the band structure parameters for strained 
In^-Gai-^As compounds as a function of alloy composi¬ 
tion, x. These ternary alloys have been modeled using 32- 
atom special quasi-randonr structures (SQSsjr^r— and 
216-atom cubic supercell (SC) structures with In and Ga 


atoms randomly distributed over the lattice sites. The 
binary alloys (GaAs and InAs) are first treated with a 
primitive 2-atom cell for highly converged reference cal¬ 
culations. The supercell approach is then adopted for 
correlated calculations of the physical properties of the 
binary InAs and GaAs, and ternary In^Gai-^As alloys 
with x = 0.25, 0.50, and 0.75. Using density functional 
theory (DFT) approach, we have obtained the composi¬ 
tional dependence for the hydrostatic deformation po¬ 
tentials of the band gap and higher conduction band 
energies at the L and X valley minima as well as two 
shear deformation potentials for the valence band split¬ 
ting at the T point. These calculations have been done 
for several random configurations of 216-atom supercell 
alloy structures to compute configurational averages of 
the physical quantities. We interpolate the DFT-derived 
compositional dependence of a physical quantity, Q, as 

Q(x) = xQ(l) + (1 - x) Q(0) - b w x (1 - x), (1) 

where 6 W is a bowing parameter, and <2(0), Q(x) and 
<2(1) correspond to GaAs, In^Gai-^-As and InAs, respec¬ 
tively. The bowing is assumed to be upward (downward) 
for & w < 0 (& w > 0). Our results predict a sizable com¬ 
positional bowing of band energies and their deformation 
potentials (6 W ^ 0), contrary to the widely used Vegard 
law (& w = 0)£ii£ The bowing effect of the deformation 
potentials is attributed to a stronger contribution of the 
p-orbitals of In cations to the formation of the alloy va¬ 
lence band states compared to that of Ga for x ~ 0.5. 

We perform the DFT calculations within the frame¬ 
work of the plane-wave pseudopotential approach com¬ 
bined with the projector augmented wave (PAW) formal¬ 
ism^ as implemented in the VASP code. 17,18 Local den¬ 
sity approximation (LDA)i^ and hybrid functional ap¬ 
proach (HSE06)2S are adopted for the DFT calculations. 
The Brillouin zone (BZ) is sampled with a 12x12x12, 
4x4x2, and 2x2x2 k-point grid for the 2-atom primitive 
cell, 32-atom SQS, and 216-atom cubic supercell struc¬ 
tures, respectively. The k-meshes are chosen to include 
the r, X and L points in the BZ. A kinetic energy cut-off 
of 400 eV is used for total energy calculations to converge 
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the total energy to 10~ 6 eV and interatomic forces to 
~ 10 -3 eV/A at least. Our supercell calculations suggest 
that Vegard’s law is a good approximation for the lattice 
constant of the fully optimized In^Gai-^As structures, 
i.e., a(x) w xai n A s + (1 - z)aG a As, where ai n A s =6.104 
A (6.032 A) and ciGaAs=5.677 A (5.611 A) are obtained 
with the HSE (LDA) functional. For the sake of compar¬ 
ison, we have also done alloy calculations using a mod¬ 
ern version of the virtual crystal approximation (VCA) 
based on the DFT approach^ This DFT-VCA scheme 
predicts, however, an anomalous compositional bowing 
of the lattice constant and band energies in contradiction 
with experiment and superccll calculations «21 Additional 
computational details can be found in the supplementary 
material (SM)<^ 

Using the lattice parameters, a(x), we have calculated 
the electronic structure of the InGaAs ternary alloys to 
obtain the relevant band energies at the T, L and X val¬ 
ley minima as a function of alloy composition as shown 
in Fig. |H The figure reveals a non-negligible downward 
bowing for the band-gap energy (E g = Er), suggesting 
that the Vegard law is a rough approximation for the 
compositional dependence of the InGaAs band gap. The 
hybrid functional calculations give a bowing parameter 
of ~ 0.55 eV in a good agreement with that of 0.43 and 
0.48 eV found in experiment at T = 0 and 300 K, re¬ 
spectively. The overall accuracy of the calculated band 
gap compared to the experimental data is in the range 
of 4-14% depending on the alloy composition and tem¬ 
perature, and whether the spin-orbit interaction is taken 
into account or not. Hereafter, we neglect the spin-orbit 
interaction, assuming that it does not significantly af¬ 
fect the compositional bowing of band energies and their 
deformation potentials. 23 This is a plausible assumption 
since the spin-orbit splitting, A so , of top valence bands 
at the F-point depends weakly on the InGaAs alloy com¬ 
position as given by A so (x), which is 0.36 and 0.38 eV 
for GaAs (x = 0) and InAs (x = 1), respectively. 

To verify that the compositional bowing does not de¬ 
pend on a particular choice of the alloy structural model 
that accounts for configurational disorder, we explicitly 
show in Fig.[T]that the 32-atom SQS and 216-atom cubic 
supercell models give consistently similar values for the 
band gap. The band gap calculated for the 32-atom SQSs 
exhibits a certain alloy broadening for the In^Gai-^As 
alloy with x = 0.5 because the three top valence bands 
are split by alloy disorder. Increasing the supercell size to 
216 atoms reduces this effect as seen in Fig. [I] This sug¬ 
gests that some physical quantities of the ternary alloys 
may suffer from a finite-size effect as a result of treating 
alloy disorder with relatively-small model structures. It 
also means that actual alloy broadening of the band gap 
is likely to be rather small in the limit of an infinitely 
(or sufficiently) large ternary alloy structure, but it may 
manifest itself in alloy nanostructures. 

Figured] also shows how the conduction band energies 
of the InGaAs alloys at the L and X valley minima de¬ 
pend on alloy composition. We find that the conduction 


Black circles: HSE, 216-atom SC 
Gray circles: HSE, 32-atom SQS 
Solid line: 6w,hse ~ 0.55 eV 
Dashed: 0 K, bw, exp ~0.43 eV 
DotDashed: 300 K, bw,exp~0.48 eV 




FIG. 1. Compositional dependence of the conduction band 
energies Er, El and Ex averaged over the one T, four L, and 
three X valley energies, respectively. The HSE-calculated en¬ 
ergies (black and gray filled circles) are defined with respect 
to the valence band energy averaged over top three valence 
bands (split by alloy disorder) at the T point. The experi¬ 
mental data are taken from Refs. [f| and [24l . Filled squares 
correspond to i?x,GaAs, which is the conduction band energy 
of GaAs exactly at the X point of the BZ. All the lines are 
given by Eq. ©■ The ’’error” bars are given by a disorder- 
induced energy splitting at the effective L (X) point of the 
virtual fee crystal BZ of InGaAs. Note that these bars are not 
related to an error, but rather quantify a real effect of alloy 
broadening. The latter are not a measure for an uncertainty 
due to different atom arrangements, which is within RT of 25 
meV only. 
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band energy, El, at the L-point exhibits strong compo¬ 
sitional bowing, suggesting that Vegard’s law is not ap¬ 
plicable for this band parameter. The conduction band 
energy, Ex, at the X valley minimum has a nonlinear 
behavior with respect to alloy composition as well, see 
Fig. [□ The figure suggests that the alloy broadening of 
these energies is more significant compared to that of the 
band gap. Unfortunately, there are no reliable experi¬ 
mental data for the higher conduction band energies of 
the In^Gai-^As alloys with x > 0 so that we provide a 
parameter-free prediction of their values. 

The higher energies are more prone to alloy disorder 
due to a different character of the electronic states at the 
L and A' valley minima compared to that of the electronic 
state at the T point. The latter is an isotropic, s-like 
state that should be rather insensitive to some structural 
anisotropy induced by disorder in the InGaAs alloy. The 
electronic states at the four L and three X valleys of the 
BZ are anisotropic due to their mixed sp-like character 
that makes them more sensitive to alloy disorder. In the 
supercell calculations, we observe it as a four (three)-fold 
degeneracy lift, i.e., energy level splitting, at the effec¬ 
tive L(X) valley in the virtual face-centered cubic (fee) 
crystal BZ, which can be defined for the InGaAs alloys 
in virtual crystal approximation. We compute the actual 
conduction band energies shown in Fig. Q] as an average 
of the corresponding four (three) energies at the effective 
L (A) valley. The alloy broadening of the band energies 
is proportional to the disorder-induced energy level split¬ 
ting, see Fig. [I] Though this splitting can be sizeable, 
the averaged conduction band energies are insensitive to 
configurational disorder with some exception for x ~ 0.5. 
That is demonstrated by using two different structural 
models of random alloys (the SQS and cubic supercells) 
that give similar averaged energies as shown in Fig. [1] 
Thus, we provide a reliable prediction for the average 
band structure of InGaAs even though alloy broadening 
might be overestimated because of the finite-size effect. 

We now describe how the conduction band energies are 
affected by hydrostatic strain. The corresponding defor¬ 
mation potential, ac, is given by Ec(e) = Eq( 0) + 3 ac e 
with C = r, L and X, where e = e xx = e yy = e zz ^ 0 are 
the only non-zero elements of the strain tensor^^ A 
larger absolute value of the deformation potential means 
a stronger strain influence on the alloy band structure, 
and the sign of ac suggests a qualitative behavior of the 
band energies with respect to the compressive (e < 0) 
and tensile (e > 0) strain conditions. Using the LDA 
and HSE functionals, we have calculated the band-gap 
hydrostatic deformation potential, ar, for the 216-atom 
cubic supercell structures of the In^Gai-^As ternary al¬ 
loys as a function of alloy composition as given in Fig. [2] 
The LDA (HSE) calculations agree (very well) with the 
experimental data for the GaAs and InAs binary alloys. 
For the InGaAs ternary alloys, a non-negligible upward 
bowing is predicted for ar, see Fig. [2] The 32-atom SQS 
model shows a similar upward behavior for the composi¬ 
tional bowing of ar, proving that we have done configu- 
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FIG. 2. Compositional dependence of the hydrostatic defor¬ 
mation potentials (ar, ar and ax) of the conduction band en¬ 
ergies that are defined as ac = [ac(e = 0+) + ac(e = 0—)]/2 
with C = F, L and X, where ac(e = 0 +) and ac(e = 0 —) are 
calculated with the HSE (filled circles and squares) and LDA 
(filled gray rhombs) functionals for the tensile and compres¬ 
sive strain conditions, respectively. Filled squares correspond 
to the deformation potential, a^, for the conduction band en¬ 
ergy of GaAs exactly at the X point of the BZ. The ’’error” 
bars correspond to alloy disorder-induced anisotropy quanti¬ 
fied as |ac(e = 0 +) — ac(e = 0 — )|. The open squares are 
related to the experimental data.— All the lines are given by 
Eq. ©. 


rational averaging of this physical parameter properly. 

Having a good agreement with experiment for the two 
binary alloys (GaAs and InAs) suggests that the HSE 
calculations should allow for reliable prediction of hy- 
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drostatic deformation potentials for the corresponding 
ternary alloys. Though three experimental data points 
for the In x Gai_ x As alloys with x = 0.0, 0.5, and 1.0 
(Fig. [2j seem to suggest a downward compositional bow¬ 
ing of the band-gap hydrostatic deformation potential, 
we would like to point out that these data are taken 
from three uncorrelated sources^ Another explanation 
is that determining a g from row experimental data re¬ 
quires separating the hydrostatic and shear strain effects 
on the band gap. This procedure relies on the knowl¬ 
edge of elastic constants for InGaAs, which are usually 
estimated from Vegard’s law. However, we have found 
that the elastic constant compositional dependence of 
the InGaAs compounds deviates from the Vegard law 
as discussed in SM , 22 and this might be a reason for the 
discrepancy. 

Besides the band-gap energy of the InGaAs com¬ 
pounds, the hydrostatic strain affects the higher ener¬ 
gies at the L and X valley minima. The corresponding 
deformation potentials are given as a function of alloy 
composition in Fig. [2] One can clearly see that there 
exists a strong upward (downward) bowing of the hydro¬ 
static deformation potential, ol (ax), for the conduction 
band energy exactly at the L (X) point. For the A' val¬ 
ley minimum, ax(% ) exhibits a non-monotonic behavior, 
reaching its maximum at x < 0.25. We find that the de¬ 
formation potential, a^, exhibits a modest alloy-induced 
anisotropy (proportional to alloy broadening) given with 
’’error” bars in Fig. [2] That is confirmed by the HSE cal¬ 
culations done for the 32-atom SQS and 216-atom cubic 
supercell structures that represent two different models 
for configurational disorder in the InGaAs ternary alloys. 
On the contrary, the hydrostatic deformation potential, 
ax, suffers from a significantly larger alloy broadening 
than a£. That is mainly due to strain-induced chemi¬ 
cal hybridization of the electronic states at the effective 
A-point of the supercell (unfolded) BZ. 

The band structure of the In^Gai-^As ternary al¬ 
loy can also be modified by two shear deformations (I) 
e = c X x e yy , Cyy = f-xx 7 ^ 0 and c zz = 0, and (II) 
e = e xy = e yx ^ 0 i 12 ’ 23 One of three top valence bands 
at the r point is not altered by shear strain whereas 
two other valence bands are split symmetrically by the 
strain. The corresponding splitting energy is given as 
A E v = D |e| in the linear approximation, which we find 
to be valid up to e ~ 0.01 at least. We define two shear 
deformation potentials as D\ = 3 b and Du = 2 x 3 1 / 2 d, 
where b and d correspond to the shear strain I and II 
according to k-p model conventionsThe microscopic 
parameters obtained from our DFT calculations can then 
be directly employed for this phenomenological model. 
The shear strain affects the band gap by decreasing it to 
E g ( 0) — | A.E7 V | /2 unlike the hydrostatic strain that can 
enlarge or narrow the gap, depending on the sign of e. 
The electronic states at the F (A', L) valley are unaf¬ 
fected (somewhat affected) by shear strain, see 

Figure 0 shows how the shear deformation potentials 
(b and d) depend on alloy composition. One can see that 
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FIG. 3. Compositional dependence of the shear deformation 
potentials, D = b or d, defined as D = [D(e = 0+) + D(t = 
0—)]/2, where D(t = 0+) and D(e = 0—) are calculated for 
the tensile and compressive strain conditions, using the LDA 
(gray filled rhombs) and HSE (black filled circles) functionals. 
All the lines are given by Eq. JT]). The recommended values 
from Ref. @ are given with arrows. The ’’error” bars are re¬ 
lated to disorder-induced anisotropy, | D(t = 0+) — D(t = 
0 -)|. 


the difference between the LDA and HSE-calculated b 
(i d) deformation potentials is of ~3% (5%) only. This 
is likely due the fact that the valence bands are usually 
well described by the LDA approximation^ In Fig. [3l 
we give some values of b (d) recommended in Ref. [6| for 
the binary alloys. Our HSE calculations agree well with 
these recommended data, predicting a somewhat higher 
(lower) shear deformation potential for InAs (GaAs) by 
~10% (0%) and ~5% (8%) for the b and d deforma¬ 
tion potentials, respectively. However, we would like to 
emphasize that there hardly exist reliable experimental 
data for a proper comparison. The present first-principles 
study does not rely on any adjustable parameters. It can 
serve as a reference for future experiments, and can be 
exploited to test empirical band structure calculations. 

Using the first-principles calculations, we can now clar¬ 
ify whether the compositional bowing of shear deforma¬ 
tion potentials is strong or not. Figure [3] clearly shows 
a nonlinear compositional dependence of the b and d 
parameters. This bowing is not as strong as that pre- 
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dieted for group-III nitride alloys in Ref. EH but it re¬ 
veals the apparent limitation of Vegard’s law, which is 
often assumed for the band energy deformation poten¬ 
tials of ternary alloys*^; We find that the composi¬ 
tional bowing of the deformation potentials is upward 
in the In^Gai-^As alloys, meaning that the strain de¬ 
pendence for the valence band splitting of ternary al¬ 
loys with x ~ 0.5 is more similar to that of InAs (not 
GaAs). A plausible explanation for this behavior is that 
the top valence bands, which show a p-like character, are 
preferably formed by the p-orbitals of In (not Ga) cation 
atoms in the In^Gai-^As alloys with x > 0.5. This is 
confirmed by comparing the relative contribution of In 
and Ga atomic orbitals to the density of states for the 
top valence bands. Finally, we notice that there exists 
some broadening of the b and d parameters because of 
the disorder-induced anisotropy with respect to the strain 
parameter sign, see Fig. [3] This broadening of d ( b ) is, 
however, weaker than (comparable to) the compositional 
bowing of d ( b). The anisotropy becomes unphysically 
large for the shear deformation potentials obtained with 
the 32-atom SQS structures (not reported in Fig. [3]) be¬ 
cause of a finite-size effect. 

In conclusion, we studied the structural, electronic and 
mechanical properties of strained In^Gai-^As ternary al¬ 


loys in a systematic manner, using a first-principles ap¬ 
proach for various alloy compositions and strain config¬ 
urations. The compositional and strain dependencies of 
the band structure parameters as well as other relevant 
material parameters (tabulated in SM, see Ref. HH) were 
obtained for these semiconductor compounds, allowing 
for device simulations within the framework of the tight- 
binding and drift-diffusion models. We found that the 
band energies and their deformation potentials can ex¬ 
hibit a significant compositional bowing. Using the en¬ 
ergy and atom-resolved density of states, the role of In 
cation atoms was revealed in relation to the bowing di¬ 
rection of the InGaAs deformation potentials. That sug¬ 
gested a simple way to understand the bowing effect in 
semiconductor alloys in a qualitative manner. 
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SUPPLEMENTARY MATERIAL 

Hybrid functional calculations 

To calculate the electronic, structural and mechanical 
properties of the strained InGaAs compounds, we em¬ 
ploy the conventional HSE06 hybrid functional func¬ 
tional with 25% of the short-range Fock exchange and 
the range-separation parameter of /i = 0.207 A -1 for all 
alloy compositions. However, ion relaxation is first done 



6 


using the semi-local GGA-PBE functional to have a good 
initial guess for atomic positions in alloy structures. Sub¬ 
sequently, we perform total energy calculations with the 
HSE06 functional to achieve a better convergence for the 
atomic and electronic structure of alloys. Even though 
this second ion optimization procedure usually affects 
the HSE-calculated electronic structure marginally, it is 
particularly important for an accurate calculation of the 
strain II deformation potential, d. The shear strain II in¬ 
duces an internal strain in the InGaAs alloy structures, 
giving rise to a piezoelectric effect. Accounting for this 
effect properly requires calculating the interatomic forces 
and band energies in a consistent manner, i.e., with the 
same functional for all the physical quantities. 

For a 216-atom supercell, we have found that the band 
structure energies are sufficiently accurate if a lower ki¬ 
netic energy cut-off of 300 eV and a T-point for the 
short-range Fock exchange part of a hybrid functional are 
adopted. That allows reducing the computational burden 
of hybrid functional calculations without compromising 
the required accuracy. The induced error (upward shift) 
is of ~20 meV, which is given by the difference between 
the 32-atom SQS (2-atom primitive cell) and 216-atom 
supercell-calculated band parameters for the binary al¬ 
loys as seen in Fig. 1 in the main text and Table SI in 
supplementary material. This shift is within an accept¬ 
able error bar since it is much smaller than the difference 
between the calculated and experimental values of the 
band energies. 

Lattice parameter definition for ternary alloys 

To understand the structural properties of the InGaAs 
compounds, we have adopted two structural models (the 
32-atom SQS and 216-atom cubic supercells) as discussed 
in the manuscript. In both cases, an effective 2-atom cu¬ 
bic cell is assumed for calculating a lattice parameter, 
a , of ternary alloys, which is derived from a simple re¬ 
lation between the supercell volume (f2) and lattice con¬ 
stant (a), f2 = N u a 3 , where N u is the number of 2-atom 
cells constituting the SQS or cubic supercell structures. 
To calculate the lattice parameter, a full structural opti¬ 
mization has been done for the ternary alloy model struc¬ 
tures. We have found essentially no compositional bow¬ 
ing of the lattice constant for the fully optimized InGaAs 
structures, i.e., 6 W ~ 0. Doing ion relaxation is crucial to 
account for the effect of alloy disorder in ternary alloys 
properly. Otherwise, it would lead to the wrong predic¬ 
tion of an upward (6 W < 0) compositional bowing for the 
lattice constant of the InGaAs ternary alloys in spirit of 
the VC A that does not allow any ion relaxation. Note 
that the HSE-calculated lattice constant for GaAs (InAs) 
is just 0.5% (0.8%) larger than the corresponding exper¬ 
imental parameter as seen in Table SI. Thus, the hybrid 
functional approach can be expected to give the lattice 
parameters of the InGaAs alloys with a similar accuracy. 

Compositional dependence of elastic constants 


Another explanation for the apparent discrepancy be¬ 
tween the DFT-calculated and experimental values of the 
band-gap deformation potential is that the measurement 
for the In^.Gai-^As alloy with x = 0.53 was done un¬ 
der general strain conditions. It means that it requires 
separating the hydrostatic and shear strain effects on the 
band-gap energy to determine the actual hydrostatic de¬ 
formation potential. This separation procedure relies on 
the knowledge of other alloy material parameters such as 
elastic constants, which are usually assumed to be given 
by Vegard’s law. However, we have found that the elastic 
constants of the InGaAs compounds exhibit some down¬ 
ward compositional bowing that might need to be taken 
into consideration for more accurate interpretation of the 
experimental data. In Table SI, we report the elastic con¬ 
stants of the binary alloys (GaAs and InAs) as well as the 
corresponding compositional bowing parameters for the 
InGaAs ternary alloy calculated with the HSE functional. 
The table shows that the elastic constants of the binary 
alloys agree very well (within ~l-2%) with the existing 
experimental data. A somewhat larger deviation of ^6% 
is found for Ci2,GaAs- 

Conduction band deformation potentials 

We would like to notice that the conduction band hy¬ 
drostatic deformation potentials, ar, ol and ax, in Ta¬ 
ble SI can be expressed as ar = r — ay and ac = 
S d , c + S„ jC /3 - a v , where C = L or X, and H^r, S d>c 
and H U| c are the deformation potentials as defined by 
Herring and Vogt [Phys. Rev. 101 , 944 (1956)], and ay 
is the hydrostatic deformation potential for the average 
of three top valence bands at the T point. The electronic 
states at the T and L (T and X) valleys of the face- 
centered cubic BZ are not affected by the shear strain I 
(II), being protected by symmetry, i.e., S u .r = 0. The 
three X (four L) valleys are split by the shear strain I (II) 
into three singlets (doublets) with a shift of A Ex = 0 and 
A E x = ±E u ,x|e| (A E L = ±2E U>L |e|/3) with respect to 
the averaged conduction band energy. The conduction 
band deformation potentials are tabulated in Table SII. 

VCA vs. supercell calculations of ternary alloys 

For the sake of comparison, we have done explicit vir¬ 
tual crystal approximation (VCA) calculations based on 
density functional theory (DFT), using a version of DFT- 
VCA proposed by Bellaiche and Vanderbilt [Phys. Rev. 
B 61 , 7877 (2000)]. In Fig. SI, the obtained results 
show that this DFT-VC A scheme predicts an anomalous 
downward compositional bowing of the lattice constant 
for InGaAs ternary alloys in consistency with findings 
of Sokeland and co-workers [Phys. Rev. B 68, 075203 
(2003)] for InGaN ternary alloys. It means that the DFT- 
VCA fails to reproduce the linear compositional depen¬ 
dence (no bowing) of the lattice parameter of InGaAs 
suggested by experiment and our supercell calculations. 
Using the DFT-VCA also results in an unphysical down¬ 
ward compositional bowing of the band gap and higher 
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conduction band energies as seen in Fig. SI. For ex¬ 
ample, the compositional dependence of the band gap 
energy has a significant contribution of the third-order 
correction (a; 3 ) in contradiction with the supercell calcu¬ 
lations and experimental data. If we impose no bowing 
for the lattice parameter in the DFT-VCA calculations 
(strained VCA), the compositional bowing of the band 


gap increases significantly, showing even stronger devia¬ 
tion from the experimental findings and supercell calcu¬ 
lations, see Fig. SI. Sokeland and co-workers have pro¬ 
posed that these failures of DFT-VCA are due to rather 
different lattice constants of alloy components (InN and 
GaN in their case, and InAs and GaAs in our case) as well 
as a very different localization of the In and Ga electronic 
states that are to be mixed in the VCA calculations. 
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TABLE I. Material parameters of the InGaAs compounds calculated using the HSE06 functional. The conduction band energies 
of the binary alloys (GaAs and InAs), Er, El, Ex, and Ex (in eV), as well as the hydrostatic and shear deformation potentials, 
ar, cll, ax, a* x , b and d (in eV), are calculated for a 2-atom primitive cell without the spin-orbit coupling included. The spin- 
orbit splitting, A so , given in the table allows accounting for the spin-orbit interaction effect on the band energies within a 
k ■ p model. In the parenthesis, we report the conduction band energies obtained with the spin-orbit coupling included; the 
experimental lattice constants have been adopted in this case. The lattice parameters, a, and elastic constants, C n, C 12 and 
C44, are given in units of A and GPa, respectively. Any material parameter, Q, for a given ternary alloy composition, x, can 
be obtained using the interpolation formula Eq. (1) in the main text, where the bowing parameters, 6 W , are obtained using 
216-atom supercell model structures for the InGaAs ternary alloys. The compositional dependence of the band energy, Ex, at 
the X valley minimum can be approximated as Ex(x) = xEx, inAs + (1 — x) Ex , GaAs — b„x(l — x) for x > 0.25. The same 
approximation can be adopted for ax- As an example, we give the parameters for the most relevant alloy, Ino. 53 Gao. 47 As. The 
experimental data are given at room temperature, Vurgaftmana et al [Appl. Phys. Rev. 89, 5815 (2001)] and K.-H. Goetz et 
al [J. Appl. Phys. 54, 4543 (1983)]. 



E r 

E l 

E x 

Ex 

A so 

ar 

Ol 

ax 

ax 

b 

d 

Cn 

C 12 

C 44 

a 

GaAs 

1.29 

1.72 

2.05 

2.09 

0.36 

- 8.54 

-3.35 

1.39 

1.76 

-2.01 

- 4.33 

117.5 

50.0 

60.3 

5.677 


(1.32) 

(1.65) 

(1.91) 

(1.95) 

(0.36) 











Exp 

1.42 

1.81 

1.97 

- 

0.34 

- 8.33 

- 

- 

- 

-2.0 

- 4.8 

119.0 

53.4 

59.6 

5.653 

Ino. 53 Gao. 47 As 

0.66 

1.57 

2.11 

2.11 

0.37 

-6.81 

-2.79 

1.50 

1.50 

-1.79 

-3.90 

96.5 

46.9 

45.2 

5.891 


(0.69) 

(1.51) 

(1.95) 

(1.95) 

(0.37) 











Exp 

0.75 

- 

- 

- 

- 

-7.79 

- 

- 

- 

- 

- 

- 

- 

- 

5.868 

InAs 

0.36 

1.58 

2.17 

2.17 

0.38 

-5.97 

-2.69 

1.45 

1.45 

-1.69 

-3.72 

85.2 

45.2 

40.3 

6.104 


(0.39) 

(1.53) 

(2.00) 

(2.00) 

(0.38) 











Exp 

0.35 

- 

- 

- 

0.39 

-6.08 

- 

- 

- 

-1.8 

-3.6 

83.4 

45.4 

39.5 

6.058 

6 W 

0.55 

0.31 

- 

0.105 

0 

-1.47 

-0.83 

- 

0.40 

-0.20 

-0.43 

15.5 

2.2 

17.9 

-0 

Exp 

0.48 

- 

- 

- 

- 

2.34 

- 

- 

- 

- 

- 

- 

- 

- 

-0 


TABLE II. Hydrostatic and shear deformation potentials of the InGaAs compounds calculated using the HSE06 functional 
for the conduction band valleys, Sd,r — ay, E d,x — ay, E U: x, Ed,z, — ay and E u ,l in units of eV, where ay is the hydrostatic 
deformation potential for the average of three top valence bands at the P point. We note that S Ui r = 0. The experimental 
data are taken from Landolt-Bomstein, Numerical Data and Functional Relationships in Science and Technology (Springer, 
New York, 1982), Group III, Vol. 17a-b. In the parenthesis and last column, we give the theoretical values of deformation 
potentials, S[) DA and ay DA , which were obtained using the DFT-LDA calculations for the binary alloys (GaAs and InAs) 
by Chris Van de Walle [Phys. Rev. B 39, 1871 (1989)], and Vegard’s law is used to estimate the value of ay for the 
Ino. 53 Gao. 47 As alloy. The compositional dependence of the deformation potential, S Ui x, can be approximated as S Ui x(*) = 
X E„, A',InAs + (!-*) H^A.GaAs ~ x(l - x) for X > 0.25. 



H<j, r — ay 

Hd,x — ay 


“d ,x ~ ay 

1—1* 

^u,X 

Sd,L — ay 

^U,L 

1—iexp 

1 'u,L 

LDA 

ay 

GaAs 

-8.54 

-1.35 

8.22 (8.61) 

-0.34 

6.30 

-8.38 

15.1 (14.26) 

19.6 

1.16 

Ino. 53 Gao. 47 As 

-6.81 

- 0.22 

5.15 

- 0.22 

5.15 

-7.09 

12.9 

- 

1.08 

InAs 

-5.97 

-0.16 

4.83 (4.50) 

-0.16 

4.83 

- 6.86 

12.5 (11.35) 

- 

1.00 

b-w 

-1.47 

- 

- 

- 0.10 

1.50 

-1.96 

3.38 

- 

- 
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FIG. 4. Compositional dependence of the lattice constant, 
a (top), band gap energy (middle), and L -valley conduction 
band energy (bottom) calculated using the HSE functional for 
two alloy models based on (i) 216-atom supercell (SC), and 
(ii) 2-atom virtual crystal approximation (VCA) by Bellaiche 
and Vanderbilt. Filled (open) black circles correspond to su¬ 
percell (VCA) calculations. Filled black triangles correspond 
to strained VCA calculations where the supercell-calculated 
lattice constant is adopted, i.e., assuming no bowing for the 
lattice constant. The solid and dotdashed black lines are given 
by the formula Eq. (1), and the dashed and dotted black lines 
are given by a modified formula Eq. (1) where a third-order 
term, — c w x 2 (l — x), is added for having a better fit of VCA 
data points. In parenthesis, we give bowing parameters for 
strained VCA. 












